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Mutual  Coupling  Effects  on  a Circular  Array  of  Dipoles 


1.  INTRODUCTION 

Studies  of  coupled  antennas  for  circular  arrays  may  be  separated  into  two  groups: 
those  that  assume  a convenient  distribution  of  current  along  each  identical  element 
regardless  of  their  relative  locations  in  the  array,  and  those  that  attempt  to  obtain 
the  actual  currents  in  the  several  elements.  Nearly  all  of  the  early  and  most  of  the 
more  recent  analyses  are  in  the  first  group,  in  which  both  field  patterns  and  impe- 
dances have  been  determined  for  elements  with  assumed  currents.  An  analysis  of 
coupled  antennas  from  the  point  of  view  of  finding  the  actual  distributions  of  currents 

for  the  N-element  circular  array  was  done  by  King1  in  1950,  and  a general  analysis 

2 

of  arrays  of  coupled  antennas  was  also  done  by  King.  This  is  done  by  formulation 
of  integral  equations  for  the  currents  that  must  be  satisfied  in  order  to  meet  the 
boundary  conditions  along  the  surfaces  of  the  elements.  The  drawback  of  this 
method  is  that  the  rigorous  solution  of  the  simultaneous  integral  equations  for  the 
distributions  of  current  in  the  elements  of  the  array  of  parallel  elements  is  very 
complicated  and  no  simple  and  practically  useful  set  of  formulas  is  obtained. 


(Received  for  publication  17  April  1979) 

1.  King.  R.  (1950)  Theory  of  N-coupled  parallel  antennas,  J.  Appl.  Phys.  21:94. 

2.  King,  R.  W.  P.  (1958)  Theory  of  Linear  Antennas,  Chapter15,  Harvard  Univer- 

sity, Cambridge, 
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In  this  report  we  shall  use  the  techniques  of  King  for  the  derivation  of  the  inte- 
gral equations  for  the  currents,  and  from  this  starting  point  use  the  computer  to 
numerically  solve  the  M nonhomogeneous  integral  equations  (using  numerical 
approximations  for  the  integrals  and  complex  matrix  algebra  techniques  to  solve 
the  simultaneous  equations).  Such  an  analysis  will  display  the  mutual  and  self 
admittance  characteristics  of  a circular  array.  We  shall  investigate  the  effects  of 
mutual  coupling  for  a circular  array  with  a single  active  element  by:  (1)  varying 
the  distance  of  the  circumferential  spacing  between  adjacent  elements  while  keeping 
the  frequency  and  the  number  of  elements  in  the  array  constant;  (2)  varying  the 
frequency  while  keeping  the  radius  and  the  number  of  elements  in  the  circular  array 
constant;  (3)  varying  the  number  of  elements  in  the  circular  array  while  keeping  the 
radius  and  frequency  constant;  and  (4)  generating  the  radiation  pattern  for  a single 
active  element  of  a circular  array.  Such  a study  will  have  an  impact  on  the  develop- 
ment of  optimum  methods  for  determining  techniques  for  radiation  pattern  genera- 
tion with  very  low  sidelobes  and  a minimum  number  of  array  elements. 


2.  MATHEMATICAL  ANALYSIS 

The  circular  antenna  array  to  be  considered  (see  Figure  1)  consists  of  identical, 
parallel,  cylindrical  dipoles  equally  spaced  around  the  circumference  of  a circle. 
Only  center-fed  dipoles  with  a half  length  h = X/4  and  a radius  a will  be  used  as 
elements.  The  center  of  each  dipole  is  perpendicular  to  the  plane  of  the  circle. 

Each  of  the  elements  of  the  array  is  thus  in  the  same  geometrical  environment.  The 
radius  of  the  circle  is  p.  The  electrical  circumference  of  the  circle  is  then  /3p, 
where  j3  is  defined  in  terms  of  the  wavelength  A by  the  relationship  2n /A. 


Figure  1.  Perspective  View  of  Circular  Array 
Antenna 
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Let  each  dipole  of  the  array  be  driven  by  a voltage  where  i = 1,  ....  m and 
m is  the  number  of  elements  in  the  circular  array.  A physically  realizable  element 
has  a base  gap  where  voltage  is  applied.  We  will  assume  a base  gap  width  of  20 
(see  Figure  2)  for  this  analysis.  From  the  geometrical  pattern  of  Figure  1 assume: 

1.  £Ja  < < 1 (there  will  be  only  a Z-component  of  current  in 
each  antenna). 

2.  The  end  faces  of  the  antenna  are  neglected,  and  the 
current  is  taken  as  zero  at  both  ends. 

3.  If  assumed  elements  are  perfect  conductors,  the  boundary 
condition  that  must  be  satisfied  is  that  the  electric  field 
tangential  to  each  element  on  its  cylindrical  surface 
must  be  zero. 


I ,RA0IUS'O 


lz(Z) fl  I 


Figure  2.  Cylindrical  Dipole, 
Showing  the  Method  of 
Applying  the  Boundary  Condition 

E = 0 at  r = a 
z 


It  is  readily  shown  that  the  currents  on  the  circular  array  of  m elements  are 
related  to  the  m base  voltages  by  m simultaneous  integral  equations  similar  to 


3.  Tillman,  James  D.  (1966)  The  Theory  and  Design  of  Circular  Antenna  Arrays. 
Chapter  1,  The  University  of  Tennessee,  Enginering  Experiment  Station. 
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m where 


for  each  K = 1 


is  the  intrinsic  impedance  of  free  space 


is  the  current  on  the  i element 


is  the  distance  from  the  i element 
to  the  m1^  element  in  the  2 plane. 


mm 


V.  is  the  base  voltage  for  element  K 


- , where  L is  the  length  of  the  dipole.  We  are  to  determine 
Substituting  for  Iz.CZ)  in  Eq.  (1)  we  get 

-j/3  «ki  v.  v 

e d"Z  = k cos /3Z  + sin/3  |z|  . (2) 


Let  lz.  = A.  cos 
these  coefficients  A. 


to  simplify  the  right-hand  expression,  then 


Using 


-j0Rki 


= cos  0 Rki  - j sin  /3  Rfci 


we  get 

m | h 

D A,  / 

i=  1 M-h 


Now.  let 


frZ  cosflR^ 

cos  T7  — dZ  " -h  C ~ 


7i~Z  sin 

17  Id 


*Rki 


d Z 


-i2n 


Vk- 


, ^ „TrZ  cos  ^Rfci 

‘lki^  = cos  T7  Hj~ 

f m ofia*Z  sin^Rki 
f2ki'Z)  “ cos  T 


~RkT 


then  Eq.  (5)  becomes 

m j h h | 

E Ai  ii  -ijh  w^j 


(5) 

(6) 

(7) 


-i2”  v 

n k 


(8) 


We  integrate  f^(Z)  and  ^ki^^  fro111  -h  to  h by  using  cautious  Romberg  extrapola- 
tion to  yield 


r 


E A,  {Flkl  - i F2kl  } - ^ Vk 


where 


Fiki-  X 'iki®dZ 

-n 


(9) 


(10) 


F2ki  = /h  '2ki®dZ- 


(ID 


Let 


Flki  ‘ •*  F2ki  = Tki 


then 


1 
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(12) 


I 


-j2n 


Ai  Tki  * 'ir  vk  • 


Now  we  shall  consider  the  M integral  equations  by  taking  a look  at  the  equations 
yielded  for  each  base  v^nage  V^(K  = 1, . . . m)  simultaneously.  We  get 

— i 277 


A1T11  + A2T12  + 


AmTlm=  V V1 


A1T21  + A2  * 22  * Am  T2m  “ ~rf 


]2ti 


A,T  . + A„T  . + 
1 ml  2 m2 


AT  = V 

m mm  T)  m 


In  a matrix  notation  Eq.  (13)  is  equal  to 


11  12 


T 


lm 


12  A 22 


. . T 


2m 


T T 

ml  m2  " 


..  T 


A1 

"*i“ 

A2 

V2 

• 

• 

• 

. 

A 

V 

m 

m 

where 


V = jjf-L  V . 
Vk  T)  VK 


Therefore, 

f 1 = V . 
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(13) 


(14) 


1 


Given  values  for  Vj,  Vg. . . . » Vm  and  having  numerically  approximated  the  integrals 

Tki(i=  1 ;K  = 1. . . . , m),  we  compute  the  values  for  the  A^'s  by  factoring 

the  matrix  T into  the  L - U decomposition  of  a row-wise  permutation  of  T and  solv- 
ing the  systems  of  Eq.  (5).  Once  the  values  for  the  A.'s  have  been  found,  we  can 
then  determine  the  mutual  and  self  admittances  of  the  circular  array.  Equation  (1) 
gives  the  relationship  between  the  currents  and  voltages  of  the  array  elements.  For 
a m- element  array. 


Y11V1+Y12V2  + ••••  + YlkVk  + ••*+YlmVm  =I1 

Y21V1+Y22V2+  + Y2kVk  +,"+Y2mVm  sI2 


Y .V,  + Y .V.+ + Y .V.  +...+Y  V =1  , 

ml  1 m2  2 mk  k mm  m m 


where  j,  ^22*  * * * ^kk*  * * * ^mm  are  se^  admittance  of  the  respective  element# 
Yab  is  1116  mutual  admittance  between  a and  b.  VK  is  the  applied  voltage  on  the  Kth 
element,  and  Ij,  is  the  current  in  the  Kth  element. 

For  an  array  with  a single  active  element,  which  for  simplicity's  sake  we  let 
be  element  number  one,  every  VK  where  K 4 1 is  equal  to  zero.  Thus,  Eq.  (16) 
becomes 


X1  =YUV1 
X2  = Y21V1 
X3  =Y31V1 


1 = Y V 

m mlv  1 


We  let  Vj  = l and  A^  is  the  magnitude  and  phase  of  the  current  on  element  K. 
Therefore, 
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A1 

II 

►< 

A2 

= Y21 

(18) 


Am=Ymr 

Thus,  the  coefficient  Aj  is  the  self  admittance  of  the  active  element  number  one  and 
the  coefficient  A„,  where  K is  unequal  to  one,  is  the  mutual  admittance  between  the 
first  and  Kl  element. 

3.  PRESENTATION  OF  NUMERICAL  RESULTS 

Mutual  and  self  admittances  were  determined  through  use  of  computer  programs 
to  numerically  solve  the  M nonhomogeneous  simultaneous  equations  [Eq.  (13)].  Fig- 
ure 3 shows  the  magnitude  of  the  self  admittance  of  the  active  element  as  the  distance 
between  adjacent  elements  varies  up  to  5 wavelengths.  In  each  case  the  frequency 
(1.  3 GHz)  and  array  size  were  held  constant.  Figures  3a,  3b,  and  3c  show  that  the 
magnitude  of  the  self  admittance  approaches  a common  value  with  increase  in  inter - 
elemental  spacing.  This  characteristic  value  was  found  not  to  be  frequency  dependent. 
Figure  4 shows  the  magnitude  of  the  mutual  admittance  between  the  active  element 
and  an  adjacent  element  as  interelemental  spacing  varies  up  to  5 wavelengths.  Once 
again  the  frequency  (1.  3 GHz)  and  array  size  were  held  constant.  Under  these  con- 
ditions the  mutual  admittance  rapidly  decreased  as  the  distance  between  elements 
increased.  Regardless  of  the  array  size,  the  mutual  admittance  value  approached 
a common  value  as  separation  increased.  As  in  the  case  for  the  magnitude  of  the 
self  admittance,  the  characteristic  value  was  found  not  to  be  frequency  dependent. 

The  self  and  mutual  admittances  were  calculated  for  a particular  array  size, 
holding  the  spacing  between  adjacent  elements  constant  while  varying  the  frequency. 
Figures  5 and  6 show  the  magnitude  of  the  self  and  mutual  admittance  characteristics 
under  these  conditions.  For  the  self  admittance  case  for  the  active  element  (Figure  5), 
it  is  seen  that  as  the  frequency  is  increased  the  self  admittance  values  for  a 2-,  10 -, 
and  50-element  array  (Figures  5a,  5b,  and  5c)  approach  a common  value.  This 
characteristic  value  was  found  to  be  independent  of  the  spacing  between  adjacent 
elements.  Figure  6 shows  that  as  the  frequency  is  increased  the  magnitude  of  the 
mutual  admittance  decreases.  This  is  due  to  the  fact  that  as  the  frsquency  is  in- 
creased the  distance  between  adjacent  elements  in  wavelengths  is  ncreased.  Again 
the  curves  approach  a common  value  independent  of  the  array  size  and  the  spacing 
between  adjacent  elements. 
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Figure  3.  Magnitude  of  Self  Admittance 
for  Varying  of  Interelemental  Spacing 
While  Holding  Number  of  Elements  and 
Frequency  Constant 
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Figure  4.  Magnitude  of  Mutual  Admittance 
for  Varying  of  Interelemental  Spacing 
While  Holding  Number  of  Elements  and 
Frequency  Constant 
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Figure  5.  Magnitude  of  Self  Admittance 
for  Varying  Frequency  While  Holding 
Number  of  Elements  and  Interelemental 
Spacing  Constant 
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Figure  6.  Magnitude  of  Mutual  Admittance 
for  Varying  Frequency  While  Holding 
Number  of  Elements  and  Interelemental 
Spacing  Constant 


Figures  7,  8,  9,  10,  and  11  illustrate  the  magnitude  of  the  self  and  mutual 
admittances  when  the  frequency  and  radius  of  the  circular  array  are  held  constant 
while  the  number  of  elements  comprising  the  array  is  varied.  As  the  number  of 
elements  is  increased,  the  spacing  between  the  elements  becomes  less  than  half  a 
wavelength  and  the  curve  becomes  linear  at  the  right.  It  was  found  that  as  the  radius 
of  the  array  was  increased,  the  self  admittance  curve  tended  to  spread  out,  that  is, 
the  second  peak  value  was  diminished  and  the  slope  of  the  linear  section  was  reduced 
(see  Figure  9).  Similar  results  were  found  in  the  mutual  admittance  between  the 
active  element  and  an  adjacent  element  for  constant  radius  and  frequency  with  vary- 
ing number  of  elements.  With  the  increased  number  of  elements  forcing  inter- 
elemental  spacing  to  be  less  than  half  a wavelength,  the  mutual  admittance  curve 
becomes  linear  to  the  right,  as  in  Figure  10.  For  different  radius  values,  as  the 
radius  is  increased  the  slope  of  the  linesr  portion  of  the  mutual  admittance  curve 
is  reduced  (less  positive),  as  seen  in  Figure  11. 

Figure  12  shows  the  coordinate  system  used  in  the  calculation  of  the  active 
element  patterns.  The  active  element  patterns  for  a circular  array  displaying  the 
effects  of  mutual  coupling  are  shown  in  Figure  13  for  a 6-,  10-,  30-,  50-,  and 
100-element  array  at  a frequency  of  1.  3 GHz  and  interelemental  spacing  of  0.  5 X . 
With  an  increase  in  the  number  of  elements  in  the  array,  there  is  a corresponding 
increase  in  the  backlobe  and  the  number  of  sidelobes.  This  is  because  unlike  a 
planar  array  whose  elements  contribute  equally  in  a given  direction,  the  elements 
of  a circular  array  point  in  various  directions. 

4.  DISCUSSION  AND  CONCLUSION 

It  has  been  shown  that  numerical  techniques  for  approximation  of  the  solution 
of  the  M nonhomogeneous  simultaneous  equations  (Eq.  (13)]  provide  useful  analytical 
data  in  the  determination  of  characteristics  of  a circular  array.  An  understanding 
of  the  effects  of  mutual  coupling  is  necessary  in  the  synthesizing  of  an  aperture 
distribution  to  give  the  best  fit  to  a specified  radiation  pattern  and  for  the  determina- 
tion of  optimization  techniques  for  the  circular  array  antenna.  Frequency,  inter- 
elemental spacing,  and  the  number  of  elements  comprising  the  array  affect  the  degree 
of  coupling  and  should  be  considered  in  the  synthesis  of  radiation  patterns  for  a 
specified  illumination  function. 
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Figure  9.  Composite  of  Figures  8a,  8b,  and  8c 
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Figure  10.  Magnitude  of  Mutual  Admittance  for 
Varying  Number  of  Elements  While  Holding 
Radius  and  Frequency  Constant 
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Figure  10.  Magnitude  of  Mutual  Admittance  for 
Varying  Number  of  Elements  While  Holding 
Radius  and  Frequency  Constant  (Cont) 
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Figure  11.  Composite  of  Figures  10af  10b,  and  10c 
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